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Abstract. In the framework of causal perturbation theory we consider a mas- 
sive scalar field coupled to gravity. In the field theoretic approach to quantum 
gravity (QG) we start with a massless second rank tensor field. This tensor 
field is then quantized in a covariant way in Minkowski space. This article deals 
with the adiabatic limit for graviton radiative corrections in a scattering pro- 
cess of two massive scalar particles. We compute the differential cross-section 
for bremsstrahlung processes in which one of the outgoing particles emits a 
graviton of low energy, a so called soft graviton. Since the emited graviton will 
not be detected we have to integrate over all soft gravitons. 



1. Introduction 

The study of the infrared problem in quantum gravity goes mainly back to the 
work of Weinberg JllJ who has investigated the infrared behaviour of virtual and real 
soft graviton emission processes. The transition amplitude for a process in which an 
emission of a single graviton occurs is logarithmically divergent. In the sum over 
an infinite number of soft gravitons emitted he claims that infrared divergencies 
in the transition amplitude cancel to every order of perturbation theory. In his 
treatment he considers virtual and real bremsstrahlung seperately. For the virtual 
ones he obtained in the sum over an infinte number of bremsstrahlungs processes 
a transition amplitude which is proportional to some power of the lower cutoff 
parameter A. So this rate will vanish in the limit A — > 0. For an infinite number 
of real processes he obtained a transition amplitude which is proportional to the 
same negative power in A. So the cutoff disappears if virtual and real processes are 
combined. From our point of view we think that the problem can not be treated 
in this rather simple way for the following reasons. The structure of the problem 
is the same as in quantum electrodynamics (QED). Here it is well known that the 
infrared divergencies cancel in lowest order but in higher orders one has to deal 
with subdivergences in an arbitrary Feynman diagram. The question of wheather 
or not these divergencies cancel to every order of perturbation theory is a topic 
which had long been investigated, see |l2|, |5| , but is still under discussion |l(| . As 
far as we know the cancellation of these divergencies is not at all clear. So we think 
that one has to investigate these processes carefully order by order. 

In this article we will reinvestigate the infrared problem in quantum gravity in 
the framework of causal perturbation theory. Especially we consider the emission 
of a real soft graviton in two particle scattering. We use the method of adiabatic 
switching with a testfunction from the Schwartz space which cuts off the interaction 
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at large distances in spacetime. This makes all expressions well defined during the 
calculation and it turns out to be a natural infrared cutoff. We therefore avoid the 
introduction of a graviton mass. The latter, although widely used in the literature, 
is less satisfactory on physical grounds because it modifies the interaction at short 
distances, too. It is not at all clear that the different infrared regularizations give the 
same result for observable quantities. To study this question we want to calculate 
the differential cross section for bremsstrahlung from first principles. 

Causal perturbation theory goes back to ideas of Bogoliubov pj and was carefully 
developed by Epstein and Glaser in the seventies. Later Scharf || applied it 
successfully to QED. The idea is to fix the first order interaction and then to 
construct higher orders of the scattering-matrix S by induction using free fields 
only. Since the free quantum fields as basic objects are operator valued distributions 
the scattering-matrix will be constructed as an operator valued functional on some 
testfunction space. The most important ingredient for the inductive construction 
is the causality requirement, which roughly states that the scattering- matrix of 
a sum of two testfunctions factorizes if the supports of the testfunctions can be 
separated in time. As testfunction space one considers for convenience the Schwartz 
space of functions of rapid decrease, since this space is invariant under Fourier 
transformation. Then all expressions are well defined tempered distributions. In 
doing this we cut off the interaction at large but finite distances which is unphysical 
in most cases and has to be removed at the end. This is the so called adiabatic 
limit where we take the limit that the testfunction goes to a constant. In this way 
we investigate the long range behaviour of the theory. 

2. Free fields 

We consider a symmetric second rank tensor field h liu (x) 1 which is a solution 
to the wave equation, i.e. = 0. All tensor indices are raised and lowered 

with the Minkowski metric r/^ with components 7700 = +1 = — Viii * = 1> 2, 3. We 
denote by kx = k ■ x = ri^k^x" the Minkowski scalar product of two four vectors. 
Throughout the paper equal indizes in one expression are always be properly con- 
tracted with the Minkowski metric irrespective of their position. We quantize the 
graviton field h} iV (x) by imposing the Lorentz covariant commutation relations 

[h^(x), h^(y)] = -ib«^D (x - y) (2.1) 

where the tensor b a ^^ v is defined by 

ha ppv = ^ v f3v + ^ocu^u, _ v »fJ v ^ (2.2) 

and Dq(x — y) is the massless Pauli- Jordan distribution. We can write down the 
Fourier- representation of the field W v ' . It is given by 

h af3 (x) = (2tt)- 3 / 2 I f k (a al3 (k) exp(-zfcx) + a" /3 (fc) t cxp(+zfc.x)) (2.3) 
3 \/2c(fc) 

Here u>(k) = \k\ and a a @(k), a a °(k)^ are annihilation and creation operators on a 



bosonic Fock-space. From (2.1) we find that they have the following commutation 
relations 

[a ap {k),a^{ky\ =b a(3 ^5 {3) {k-k') (2.4) 

It should be emphasized that this quantization procedure is only formal in the sense 
that we don't describe physical gravitons, since our contains too many degrees 
of freedom for a real gravitational field. There must be imposed further conditions 
like a specific gauge condition, etc. in order to obtain physical gravitons. This 
was done in B. The massive scalar field (f(x) obeys the Klein-Gordan equation 
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(□ + m 2 )ip(x) = 0. It has the following representation in terms of annihilation- and 
creation-operators : 

<p{x) = (2^)- 3 / 2 f d J— (a(p) exp{-ipx) + ~a{pf exp(+^x)) (2.5) 

J y / 2uj{p) 



The energy to(p) is as usual given by uj(p) = \J p 2 +m 2 . The operator a{p) an- 
nihilates a particle with momentum p whereas a(p)' creates an antiparticle with 
momentum p. These operators obey the following commutation relations 

[a(p),a(p^]=S^(p-p') 

and all the other commutators vanish. 

3. Interaction of massive scalar matter with gravitation 
The energy-momentum tensor of </?(#) is given by 



(3.1) 

The interaction with the gravitational field will be described through the following 
first order coupling 

T 1 {x)= % -K:h^{x)b aPlxv T^{x) : (3.2) 



where k = V327rG with Newton's constant G. This leads to the following explicit 
form of the first order interaction 

Ttix) = 1 -k\2 : h a0 (xMx)lip(x),p : -m 2 : h(x)<p(x? <p{x) :] (3.3) 

where we have introduced the trace of the graviton field h(x) = h(x)^- 

4. 5-MATRIX TO THIRD ORDER 

For the calculation of the bremsstrahlung process we need to know the explicit 
form of the S*-matrix up to third order in the coupling constant k. This means 
that we have to construct the time-ordered product Ts(xi, x%, x^). According to 
the inductive construction of Epstein and Glaser we proceed by first calculating 
Tiixi, X2) which will be done in the first subsection and after that we calculate the 
relevant terms, i.e. terms with four external scalar field operators and one graviton 
field operator, in the next subsection. 

4.1. Time-ordered product to second order. For the const ruc tion of the time- 
ordered product T% (£1,2:2) we start with the first order term (|3.3[) and apply the 
inductive construction due to Epstein-Glaser as in Scharf's book |8|. First of all we 
have to build the two distributions 

A' 2 (x 1 ,x 2 ) = -T 1 (x 1 )T 1 {x 2 ) (4.1) 

and 

R 2 (x 1 ,x 2 ) = -T 1 (x 2 )T 1 (x 1 ) (4.2) 

Since A' 2 and R' 2 are products of normally ordered field operators we apply Wick's 
theorem to obtain a normally ordered expression. We are interested in scattering 
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processes so we have to concentrate on terms with zero or one contraction only. 
The contraction of the free field operators are defined by 



h af) {x 1 )h^{x 2 ) := - ib a ^ v D i ^\x 1 - x 2 ) 
ip{xi)ip{x 2 ) := - iD^\x\ - x 2 ) 



(4.3) 



where (x) is the positive frequency part of the massive Pauli- Jordan distribu- 
tion. For the splitting operation we need the difference 



D 2 (x!,x 2 ) = R! 2 (x-t,x 2 ) - A' 2 (x-t,x 2 ) 



(4.4) 



This distribution has causal support and it is given by 



D 2 {x u x 2 ) = *y [ 4 : ^(x^h^ix.Mx^^ix^l : d^d^D m ( Xl - x 2 ) 
+ 4 : h^ix^ix^ix^Mx!),* : dfdfD m (x x - x 2 ) 
- 2m 2 : h^ix^hix^x^^x^ : 3^D m {x x - x 2 ) 
h af3 (x 2 )h(x 1 Mx 2 )li P (x 1 ) : dfD m {x x - x 2 ) 



- 2m z 

- 2m 2 

- 2m 2 



h{x 2 )hP a { Xl )^{x 2 )ip{ Xl )] p : d^D m ( Xl - x 2 ) 
h{x 2 )h' m {x 1 )ip(x2yip{x 1 ) XT : d x p l D m { Xl -x 2 ) 
+ m 4 : h(x 2 )h(xi)ip(x 2 )ip(xi)^ : D m (xi - x 2 ) 
+ m 4 : h(x 2 )h(xi)ip(x 2 )^ ip(xx) : D m (xi - x 2 ) 

+ 4 : v(x 2 )l<p(x 2 ),Mxi)l<P(xih ■ b a ^D ( Xl - x 2 ) 
+ 2m 2 : ^(x^^x^^xi)^ ip(xi) : D (a:i - x 2 ) 
(p{x 2 y ip(x 2 )ip(xi)l x ip(xi)„ : D (xi -x 2 ) 

<f(x 2 y \p(x 2 )ip(x 1 )' ! ' ip(xi) : Doixx - x 2 ) 



+ 2m 2 
-Am 4 



(4.5) 



To obtain the time-ordered product T 2 we have to split the numerical distribution 
in D 2 according to it's singular order (see the Appendix for the definition of the 
singular order). The singular order of the Pauli- Jordan distribution is u> = — 2, so 
the splitting is trivial and we can use the formula 



D m {x)=D r *(x)-D%(x). 



(4.6) 



The retarded distribution R 2 {x\,x 2 ) is then given by D 2 {x\,x 2 ) if we replace all 
Pauli- Jordan distributions by the retarded ones D ret . Finally the time-ordered 
distribution T 2 (x i,x 2 ) is given by 



T 2 (x!,x 2 ) = R 2 (xi,x 2 ) - R' 2 {xi,x 2 ) 



(4.7) 
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For later reference we give here the explicit expression 

^2 



T 2 (xi,x 2 ) = i— 



4 : h^(x 2 )h^( Xl Mx 2 )M^)l ■ d^d^D F m ( Xl - x 2 ) 
+ 4 : h af) (x 2 )hP' J {x 1 )^{x 2 )}Mx 1 ), a : df^D F (xi ~ x 2 ) 

- 2m 2 

- 2m 2 

- 2m 2 

- 2m 2 



h a ^x 2 )h( Xl Mx 2 )^(xiy ■ d x a 2 D^{ Xl -x 2 ) 
h^(x 2 )h(x 1 Mx 2 )l l p(x 1 ) : dfD F ( Xl -x 2 ) 



(4.8) 



2m 2 

V 



h{x 2 )hP"(x 1 )ip{x2Mx 1 )l : - x 2 ) 

h(x 2 )h^( Xl Mx 2 )^( Xl ), a : d^D F ( Xl -x 2 ) 
+ m 4 : h(x 2 )h( X i)ip(x 2 )(p(xi) i i : D^xi - x 2 ) 
+ m 4 : h(x 2 )h(xi)ip(x 2 )^ ip{xi) : D^xi - x 2 ) 
+ 4 : ip(x 2 )l p <p(x 2 ),Mxi)l<P(xih ■ b app °D F ( Xl - x 2 ) 
+ 2m 2 : ip(x 2 )^ p (p(x 2 ), p ip(xi)^ ip(xi) : D F (xi - x 2 ) 
(fi{x 2 ) ] '(p(x 2 )ip(x 1 )'] a (p(x 1 ), a : D F (xi -x 2 ) 
</?(x 2 ) t ip(x 2 )ip(xi) J ' ip(xi) : D F (xi - x 2 ) 
-A: h^ix^ix^lvix^phP^x^ix^ix^ : 
+ 2m 2 : h af3 (x^tpixi^tfix^^hix^ipix^ tp(x 2 ) : 
+ 2m 2 : h(xt)(p{xif \p(xi)h pa (x 2 )ip(x 2 )] p ip(x 2 ), a : 
-m 4 : h(xi)ip(xiy ' ip{xi)h(x 2 )ip{x 2 y ip(x 2 ) 
This result will be needed in the next subsection for the calculation of T 3 



4.2. Time-ordered product to third order. In the previous subsecti on w e have 
calculated the distribution T 2 ( 18 ) from the given first order coupling T\ ( 3^3 ) . Now 
to obtain the time-ordered product to third order we proceed in much the same 
way as in the calculation of T 2 . We calculate the distributions R' 3 (xi, x 2 , x 3 ) and 
A' 3 {x\, x 2l x 3 ) which are given by 



R' 3 (x 1 ,x 2l x 3 ) = T 2 (xi,x 3 )Ti(x 2 ) +T 2 (x 2 ,x 3 )T 1 (xi) + Ti(x 3 )T 2 (xx, x 2 ) (4.9) 

(4.10) 

(4.11) 



and 

A' 3 {x 1 ,x 2 ,x 3 ) = Ti(x 2 )T 2 (xi,x 3 ) + T 1 (xi)T 2 (x 2 ,x 3 ) + T 2 (xi,x 2 )Ti(x 3 ) 
where T\ and T 2 arc defined by 

Tx(xi) := -T^Xi), i = l,2,3 



and 

T 2 (xi,Xj) := -T 2 (xi,Xj) + T 1 (x i )T 1 (x j ) + T 1 (x j )T 1 (x l ) 1 i,j = 1,2,3 (4.12) 

It should be noticed that we only collect those terms which have two contractions 
involving all the spacetime points x\, x 2 and x 3 , i.e. there is no product of contrac- 
tion functions with the same argument. After we have built the causal difference 
D 3 = R' 3 —A' 3 we split the numerical part into R 3 with support in the backward light 
cone and A 3 with support in the forward light cone. The time-ordered distribution 
T 3 is then given by 



T 3 (xi,x 2 ,x 3 ) = R 3 (x ll x 2 ,x 3 ) - R' 3 (x 1 ,x 2 ,x 3 ) 



(4.13) 
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Because the calculation is staightforward and not very enlightning we give only the 
result 



T 3 (x 1 ,x 2 ,x 3 ) = ^- 



8[: h^(x n ( 1 ))ip(x n ( 1 ))l r tp(x n (2))ty<p(Xn(2))eV(x*(3))<a : 



+ : /i pCT (x 7r(1) )^(x 7r(1) ), (T ^(a; ;r( 2))J 7 ^(x 7r(2) ), £ ^(x 7r( 3 ) )J Q : 
x b a & e D£(x <2) - x A3) )d;^d x p ^D^{x A1) - x <3) ) 
+ 4m 2 : h pa (x n{1) )ip(x n{1) )] p ip(x v{2) ^(p(x n{2) )(p(x n{3) ), a : 

+ : h pa ' {x^ { i))Lp{x^ { i ) ), p Lp{x^ {2 ))^{x^(2))^{x^(Z))n ■ 
x Do(x n ( 2 ) - x v (s))d^ {3) dl" (1) D^(x vW - x n{3) ) 
+ 4m 2 : /i p ' T (x 7r(1) )(^(x 7r(1) )^(x 7r( 2))^(x 7r(2) ) e (^(x 7r( 3 ) ) : 

+ : h p ^(x wW )ip(x„ w ), p cp(x„^)lcp(x % ^)^(x^ 3 ^ : 

X D^(X^ 3 ) - X^ 2 ))d^ {1) D^(x w (i) - ^(3)) 

+ 4m 2 : h(x w{1) )<p(x n ( 1) y<p(x n ( 2 ))tyV(Xn(2))>e<P(XT r (3))0 ■ 

+ : h(x vW )(p(x vW )(p(x v ^))ly^(x v (2))^(x v (3))]f3 ■ 
x b a *> e D£(x n{3) - x A2) )dl^ Di(x nW - x v{3) ) 
+ 2m 4 : h(x n{1) )ip(x^ 1) ) f ip(x v{ 2))ly<p(x w (2))^(x w ( 3 )) ■ 

+ ■ ^•(x 7r (i))(^(x 7r(1) )( / 5(x 7r(2) )t r (^(x 7r( 2))^^(x T (3)) t : 

X £>0 ( x tt(3) - X w ( 2 ))D^(x w{1) - ^(3)) 

+ 2m 4 : /i(x Tr ( 1) )(^(x 7r(1) ) t ( / 5(a; 7r (2))V(a ; x(2))¥'(a; T (3))./3 : 

+ : M^(l)M^(l)M^(2))V(a7r(2)Mz7r(3))|g : 

x D£(x M3) - x w{2 ))d x ^ (3) D^(x w{1) - x M3) ) 
+ 8m 4 : /i pff (a; 7r(1) )(^(x 7r(1) )^(a; 7r( 2 ) )V(a;7r(2))¥'(a;7r(3)) = 
+ : h pa ' (x nil) )ip(x n(1) ), p ip(x ni2) )' i \p{x A2) )ip{x A3) )^ : 

X Dq(x^ 3 ) - X w(2 ))9ff' r(1) D^(x w (i) - X 7r ( 3 )) 

+ 4m 6 : /i(x T(1 ))(^(x 7r(1) )V(x 7r (2))V(a ; 7 r(2))¥'(x 7r (3)) : 
+ : ^(x 7r (i))(y5(x 7r(1) )( y 5(x 7r(2) )V(^(2))^(x 7r (3)) t : 

x Do(x n ( 3 ) - x n{2 ))D^(x w{1) - x M3) ) 

(4.14) 

Due to the sum over all permutations of the indices this expression for T 3 is ob- 
viously symmetric in it's arguments as it is required by the definition of the time- 
ordered product. The third order S'-matrix for the brcmsstrahlungs processes is 
then given by 



d 4 xi . . . d 4 x 3 T 3 (xi,x 2 ,x 3 )g(x 1 )g(x 2 )g(x 3 ), g € 



(4.15) 
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5. TWO PARTICLE SCATTERING WITH BREMSSTRAHLUNG 



5.1. Towards the differential cross section Step 1: 5-matrix element. In 

the preceeding section we have shown explicitly the construction of the time-ordered 
product to second and to third order. The S-matrix to third order was given by 
(4.15). Now we consider the scattering process of two massive scalar particles in 
which a bremsstrahlungs graviton is emited in the final state. That is we want 
to calculate the expectation value of S 3 (g) between the following initial and final 
states 



= / dViAl^iCpO^^MpOVgi)^ (5.1) 

9 f = j d 3 ^d 3 ftd 3 S^(ft,^)^ , '(£)a'»'(£)t a (g ! )t a (^)tn (5.2) 

These are wave packets in momentum space where ipi, i = 1,2 are one particle 
wave functions from L 2 (M 3 ), ^ is a two particle wave function from L 2 (R 3 R 3 ) 
and ft*" is a tensor- valued square integrable wave function on R 3 . The vector SI is 
the vacuum vector in the asymptotic Fock-Hilbert space. We want to calculate the 
following S- matrix element: 



S fi = ^ f ,S 3 (g)^) 



(5.3) 



To show the details of the calculation we restrict ourself to the following term from 
the above computed T 3 (4.14) 



T 3 l >(xi,x 2l x 3 ) := in 3 : h pa (xi)ip{x 3 ),0ip(xi)] p ip{x 2 )l y ip{x 2 ), l 
b a ^D^(x 2 - x 3 )d^d^D^( Xl - x 3 ) 



(5.4) 



where the (1) refers to this particular term of T 3 . The ^-matrix element corre- 
sponding to this contribution is then given by 

Sf) = J d 3 p 1 d 3 q{d 3 kd 3 p 2 d 3 q 2 d 4 x 1 . . .^XsMP^Mqi^iAi&W^k) 

x 5 (x 1 )... ff (x 3 )(a^(fc) t a(p 2 ) t a(^) t O,T 3 (1) (x 1 ,a ;2 ,x 3 )a(p 1 ) t a( ( f 1 ) t r!) 

(5.5) 

Then we have to compute the vacuum expectation value 

(Sl,a{q 2 )a{f 2 )a» v {k) : h^faMx^yix^fa^fa)* : o(pl)t (gl)tn) (5.6) 

This can easily be evaluated by inserting the Fourier representations of the free field 
operators ( [2.3|) and (2.5) and the use of the commutation relations in momentum 
space (2. 4), (2. 6). The result is given by 



(SI, . . .Si) = (27r)- 1B /26^^( 3 2w(^)w(ft)w(ft)w(fi)w(&)) 1/2 

exp[+i(kx! -p\X 3 +q 2 X! +p 2 x 2 ~ qix 2 )] Piq 2 p 2 qf 

+ exp[+i(kx! - qix 3 + q 2 xi + p 2 x 2 - piX 2 )] q{ q 2 P 2 p\ 
+ exp[-H(fcri - pix 3 + p 2 x x + q 2 x 2 - qix 2 )] pj p p 2 qT, q\ 
+ exp[+i(kxi - qix 3 +p 2 xi + q 2 x 2 - pix 2 )] q{ P 2 q 2 p\ 



(5.7) 
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The ^-matrix element then reads 



Sf) = ^(27^)- 15 /Vy d 3 p 1 d 3 q 1 d i kd i p 2 d 3 q 2 d i x 1 . . . d*s 3 &(?l)^($i)*(0!, fc) 
r"(k)g( Xl ) . . . g(x 3 )b a ^D^(x 2 - x 3 )8? &? D^fa - x 3 )b^ 



We introduce the abbreviations h . . . h for the four parts of . Then we consider 
the spatial integrations in I\. They can be carried out by inserting the Fourier 
transforms of the testfunctions as well as the Feynman propagators: 



h= / d i x 1 ...d i x 3 g(x 1 )...g(x 3 )D^(x 2 -x 3 )d^d^D^ n {x 1 ~x 3 ) 



exp[+i(fca;i - p x x 3 + q 2 xi + p 2 x 2 - qix 2 )] 

= (2tt)- 10 J d 4 Xl . . . d 4 x 3 d 4 h . . . d%g(h) . . . g(l 3 ) exp[-i(l lXl + l 2 x 2 + l 3 x 3 )\ 

Dq (l 4 )ex.p[-il 4 (x 2 - x 3 )]d^d^D^(l 5 )ex.p[-il 5 ( Xl - a*)] 
cxp[+i(fca;i -piX 3 + q 2 xx + p 2 x 2 - qix 2 )] 

= (27r) -10 / d i x 1 . . .d i x 3 d i l 1 . . .d 4 l 5 g(li) . . . g(l 3 ) exp[-i(hxi + l 2 x 2 + l 3 x 3 )] 



= (2tt) 2 / d%...d%g(h)...g(l 3 )D^(l 4 )l^b^(l 5 )6^(h+l 5 -k-q 2 ) 



(32w(k)w(pi)w(p2)u)(qi)u>(q 2 )) 



exp[+i(fca;i - pix 3 + q 2 x 1 + p 2 x 2 - qix 2 )] pf q 2 P 2 Qi 
+ exp[+i(fca;i - 910:3 + q 2 xi + p 2 x 2 - pix 2 )] q{q 2 p\p\ 
+ exp[+i(fca;i - p\X 3 + p 2 x x + q 2 x 2 - qix 2 )] pf p% ql q\ 
+ exp[+i(fca;i - q x x 3 + p 2 X\ + q 2 x 2 - pix 2 )] q !i x p 2 q 2 p\ 



(5.8) 



exp[-il 4 (x 2 - x 3 )\\Hl D^(h) cxp[-i/ 5 (xi - x 3 )] 
cxp[+i(fca;i - p\x 3 + q 2 x\ + p 2 x 2 - qix 2 )] 



S {4) (l 2 + U-P2+ qi)5 {4) (h -k- h + Pi) 



(5.9) 



We can do the integrations w.r.t. l 4 and I5 immediately and get 




(5.10) 



INFRARED BEHAVIOUR OF MASSIVE SCALAR MATTER COUPLED TO GRAVITY 

Then the first term of becomes 

Sf^ = (2vr)- n / 2 z K 3 / ^dVirf 3 9ld 3 ^V2d 3 ^^i(pi)^2(<fi)*(p2,g 2 )0 All/ (fc) 
J d% . . . d%g(h) . . . g(l 3 )b al3 ^r)F (h+l 3 + Pl -k- q 2 ) 

[-1? -p1 + k a + qq] [-11 -pl + k°+ q%}Di(-h + k + q 2 ) 

6^?^?f(32o;(fc)a;(p 1 )a;(p 2 )a;(gi)a;( g l))- 1/2 

S {i) (h + l 2 + h +Pi+qi-q2-P2-k) 

d% . . . d%g(h) . . . §(!s)v#* T d±^±±f±® v»» 

{h + h+pi-k- q%y + tO 



m? — (—li + k + q 2 ) 2 — iO 1 

<5 (4) (k + h + h +Pi + qi-q2-P2-k) 

(5.11) 

where we have inserted the propagators explicitely. In this expression we can carry 
out the adiabatic limit in the variable l 2 since it doesn't appear in the argument 
of the propagators. We do this in the following way. First of all we choose a fixed 
testfunction 50 € 6>(K 4 ) with the property go(0) = 1. Then a new testfunction is 
defined by 

9e{x) := go(ex) (5.12) 

The adiabatic limit g e — > 1 is then equivalent to the limit e — ► 0. The Fourier 
transform of g £ is given by 

&(P) = ^So(f) (5-13) 

lim Igog) = {2*) 2 5^{p) (5.14) 



Then we have 



Using this in ( 5.11 ) the integration w.r.t. l 2 gives just the factor (27r) 2 . 

We are interested in the S'-matrix element in the limit where only gravitons at 
low momenta (soft gravitons) are emited. So we omit all small quantities in the 
numerators, and we take only the leading singularities in the denominators. This 
is known as the eikonal approximation in the literature, see e.g. H. With this 
approximation we obtain for the argument of the massless propagator 

(pi - 92 - k + l x + l 3 ) 2 = ( Pl - q 2 ) 2 + (-k + h+ h) 2 - 2{pi - q 2 )(k - h - l 3 ) 
= (pi - q2) 2 + o(e) 

(5.15) 

where e is the scaling parameter from the adiabatic limit. The argument of the 
massive propagator becomes 

m 2 - (-h + k + q 2 f = m 2 - <£ - (-li + kf + 2q 2 (h - k) 

(5.16) 

= 2q 2 (h-k) + o(e 2 ) 
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since the momentum q 2 is on the mass shell. Then the S'-matrix element becomes 

S { f ) Jl) = (2^)- n / 2 ZK 3 J dVld 3 ?lrf 3 fcrfV 2 rf 3 ^^l(Pl)^ 2 (gi)*(p2,^)^(fc) 

f*h*bKhmr»- <*:§ ^--lr " (..it) 

Pi &vl<A (32a;(fc)w(pi)w(p 2 )w(5i)w(q 2 )) 1/2 
<5 (4) (h + k + Pi + qi - <?2 - P2 - k) 

Now the adiabatic limit in l 3 can be done as above and we arrive at the result for 
the S'-matrix element 

fig'' - (2TT)-y 2 ZK 3 J dVld 3 9lrf 3 ^V2d 3 fe^l(A)4(ft)*(pl,fe)^(fc) 

Jdhg(li)b {pi _ q2)2 2q2{h _ k) b (5.I8) 
pf9 2 J' 2 Q , i(32w(A:)cj(pi)cj(p 2 )a)((7i)w((f 2 )) 1/2 
(5 (4) (h + pi + qi - q 2 - p 2 - k) 

There remains to calculate the other three parts S^ 1 ^ . . . S^ l4, \ The calculations 
are similiar to the one above so we only present the results 

i l5llj (qi-q2) 2 2q 2 (li-k) b 
qiq P 2Plp\{^w{k)^{pi)u{p 2 )u{qi)u{q 2 )) 1/2 
<5 (4) (h + qi + pi - q 2 - p 2 - k) 

/" f/ 4 /,fifMh a/ ^ £ ^2 "Pi) (P2-Pl) yivpo 

Jdhg(li)b {pi _ p2) 2 2p2{h _ k) b (5.19) 

PlP292 9l( 32w W W (PlV(P2)w(<?l)w(<f 2 )) _1/2 

<5 (4) (ii + qi + pi - 92 - P2 - k) 
jai l9W t> {qi - P2) 2 2 p 2 (i l -k) b 

gfp29 2 Pl( 32w (^) W (Pl) W (P2)w(<fl)w((f 2 )) V2 

<5 (4) (ii + 9i + Pi - q2 - P2 - k) 

Now we are ready to consider the differential cross section. This will be done in 
detail in the next subsection. 

5.2. Towards the differential cross section Step 2: Calculation of \Sfi\ 2 . 
We have to calculate the absolute square of the S-matrix clement sj.^ . We do this 
again separately for the four parts S^' 7 ^ . . . sj*' /4 '. For the first part we start with 
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( 5.1 8| ) and write 

|5^ Jl) | 2 = (2tt)- 7 k 6 J d 3 p 1 d 3 q 1 d 3 kd 3 p 2 d 3 q 2 d 3 p^d 3 q^d 3 k / d 3 p^d 3 q 2 ' 

/ i a f at \ r i (J i o \ 

Vxmi W 2 -p[T 2 q >(i 2 -k>) b 

v'i 12 P ' P2 1 l'i (32w(fc)w(pi)w(p2)w((fi)w(g2)) 1/2 
(32uj(k')uj(pi)uj(p 2 ')uj(q{')uj(q 2 ')) 1/2 

<5 (4) (Zi + Pi + qi - Q2 - P2 - k)6W (h + p'x + q{~q 2 ~P 2 - *') 

(5.20) 

where a * indicates complex conjugation. We use the completeness relations of the 
final states, i.e. 

X)*(P2>«5)*(P2'.<5')* = 5&(f 2 -/T2 , )5 (3) («2 - q 2 ) (5.21) 
/ 

Y^fr v {kW'"' {k'Y = rf v Tf' v 'S^\k - k') (5.22) 
/ 

where the sum is taken over all final states. Then we have 
Mp1'TM<B'T I d%d%g(h)g(l 2 )*b a ^ e 



f 



^sitf-p?) (q p 2 -p p i) 



{q 2 - Pl) 2 2q 2 (h -k) 



Pi™* 6 fe-K) 2 2q 2 (l 2 -k) Pl <Z2P2<Zl 

+ Pl + qi - q 2 - p 2 - k)6W (h + p[ + q[ -q 2 -p2-k) 

yj [32uj(k)cj(pi)uj(p2)Lj(qi)uj(q 2 )) (32uj(k)uj(p 1 ')uj(p' 2 )uj(q{')uj(q2)) 

(5.23) 

Now we assume that the wavefunctions and ip 2 are sharply concentrated around 
the initial momenta pi and Qi, so that we can replace the momenta in the propaga- 
tors by these initial momenta. Then the expression can be simplified to 

J2\ S fi' h) \ 2 = (2n)- 7 K e J d%dn 2 g(h)g{l 2 )* J d 3 kd 3 p 2 d 3 q 2 b a ^ e 

(92 ~ Pf ) (it - Pi) 3 r . P , r! e h a'0We' (l2 ~ pf) {4 ~ A) 
{q2-Pi) 2 2q 2 (h - kf l q2P2q * (q 2 -Pi) 2 2q 2 (l 2 - k) 

Pi 1 2 P 2 <lt [32uj{k)uj(pl)uj(p' 2 )uj(q l )uj(qi)] 1 



S (4) (h +Pl+qi-q 2 -p 2 - fc)5 (4) [h + p[ + q[ - q 2 - P2 - k) 

(5.24) 
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We observe that the last integral herein depends on the initial state only. We denote 
it by F(P, li, l 2 ), where we have introduced P = p 2 + q 2 + k. We have 



F{PMM) = / rfVidVirfVi'^ftViCpi)^^)^^!')*^^')* 

<5 (4) (h +Pi+qi-q2-p2- fc)5 W (h + p'i +q[-q2-p2-k) 
\-8 / j4„. .4„, 73-. j3„-'?. \„7. /„r>'\*„". /f„-*'\* 



= (2tt)- 8 / d%id*y 2 d 6 pi...d d qrMPi)Mqi)i>i(Pi TMQiT 

exp [-i(Z a + pi + <7i - P)yi] exp [+i(i 2 + pi + ?i - -P)^] 
= y d 4 yid 4 y 2 V , i(yi)'02(yi)V'i(2/2)*V , 2(2/2)*exp[-^ 1 2/ 1 ] exp[+^ 2 y 2 ] 
exp[-i(yi - y 2 )P] 

(5.25) 

The function F is normalized according to 

J F(P,l l7 l 2 )d i P^(2n) 2 J <fiy\Mv)\ 2 \Mv)\ 2 exp[-Kh-h)v] (5.26) 

and it is concentrated around P = p 2 + q 2 + k ~ pi + qi. In the limit of infintcly 
sharp wave packets we may represent it by 

F{P,hM) = (2n) 2 S^(P- Pl - qi ) J d 4 2/ ^i( 2 /)| 2 |V2(y)| 2 exp[-z(/ 1 - l 2 )y] 

(5.27) 



This will be inserted in (HH) and we obtain 



X]l4i' 7l) | 2 = j d%d%g(h)g(l 2 )* J d 5 kd?p 2 d*q 2 b a ^ 

(Q 2 -P?) (q 2 -Pr) a P Me' iQz -P?') (g| 

(92 - P*) 2 2q 2 (Z 1 - fc) 4 ^ 2 ^ (ga - p 4 ) 2 2q 2 (l 2 - k) 

Pi <1 2 P 2 4 [^2uj(k)uj(p' l )uj(p 2 )uj(ql)uj(q 2 )] 1 

5^(p 2 + q 2 +k-p i -q i ) J d 4 y|Vi(y)| 2 |^2(2/)| 2 exp[- l (? 1 -/ 2 )y] 

(5.28) 

For the definition of the cross section we follow sect. 3.4 of ||. We can set l\ = l 2 
in the last integral because we neglect all contributions o(e) in the numerator, see 
the remarks after ( |5.11 ) . With the notation as in || we define the cross section as 



a = hm irR 2 ^ \S fi {R)\ 2 (5.29) 



f 



where R is the radius of the beam of incoming particles. Then we restrict ourself 
again to (5.4) of T3 and furthermore to the first integral I\ ( |5.9D of S/j. If we denote 
the integration variables p 2 and q 2 by pf and qf, respectively, then we obtain for 
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the cross section 



-i w(pi)<j(qi) 



d 3 k d 3 q) tPpf ua0ye (q?-P?) ^-P?) 3 



d%d%g(h)g(l 2 y 



2u(k) MQf) Mff) " (<?/ ~ P,) 2 2g/(li - ff / p / ?1 

(g/ -P°) (g/ -Pi ) a' p' y £ ^ (3 ) f - - + a ) 

(q f - Pi y 2q f (l 2 - k) Pl q f P f % W + «/ + * V* *J 

+ w(g» + - w(pi) - uj(ql)) 



(5.30) 



To carry out the remainig integrations we have to know the integrand explicitely. 
So at this point of the calculation the tensor structure of the integrand becomes 
important and we have to write down the various contractions between the four- 
vectors pi,qi,pf,qf. With the tensor 6 M1/p<T ( $L2\ ) we get for the cross section 



(2tt)- 5 k 6 



d%d%g(h)g(l 2 y 



d 3 k d 3 q}- d 3 p) 



((p. ? .)2 _ 7^4) J J 2u(k) MQf) Mp~f) 

u){qj)uj{f f ) - q~fff)(piqi) + (u(dj)u(q^) - qf$)(w(fi)u(ff) - pip}) 

~ ( w (g/) w (Pi) - fffi)(u(Pf)u(Qi) -P"f0i) - {u(fi)u(ff) -Pip~f)(piqi) 

- (PiQi) (w(pi)w(p)) - pip)) + pf (uj(p})uj(ql) - p f qi) (u)(q~}) 2 - qf) 

- (w(g/)w(p» - q"fff) {p % qi) - (uj(qf)Lj(qi) - q}ql) (w(pl)w(p/) - p t pf) 
+ Mg/Mpl) - 4}Pi)(u(jPf)u(qi) - + (uj(pi)u){p f ) -p~iP~f)(piqi) 

+ {Piqi)(u{ft)uj(pf) - pip}) - p 2 {uj{p f )uj{q l ) - ffql) (u(q~f)w(fi) ~ q}pl 

5^{pjf +q~f+k-pi - ql)5(uj(p}-) + Lj(q~}) + Lo(k) - u>(pl) - u(ql)) 



(q f -pi) 2 2q f (h - k){q f - Pi f2q f {h - k) 



(5.31) 
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It is convenient to choose the center-of-mass system, defined by p\ + ql = 0. We 
can easily do the integration w.r.t. pfi 



1 (2tt)- 5 k 6 
a — 7" 



((Pi Qi 



d 4 hd%g(h)g(l 2 y 



d 3 k d 3 



If 



1 



2u>(fc) M<ff) 2u{q { + k) 



(uj(qy)uj(q) + k) + q f (q f + fc))0i?i) + (uj{q' f )uj{q 1 ) - q}ql)(uj(pl)uj(qf + k) 

+ Pi((ff + k)) - (u(q~f)uj(pi) - fffi) (w(q~f + k)uj(ql) + (q) + k)qt) 

- (uj(pl)uj(q) + k) +Pi{qj + - {p l q l )(uj{pl)uj{q} + k) + fi(q) + k)) 

+ Pi 2 (v(q~f + k)u)(q}) + (q} + k)ql) (u(q}) 2 ~ q} 2 ) ~ (u)(q})u{q~f + k) 

+ <7/(<7? + k )) (Pm) - - + k)+ Pi(0 + k)) 

+ (u(q})u(pl) - q~}p~l) (uj(q} + k)ui(qi) + {q} + k)ql) + (ui(pl)u(q'} + k) 
+ Pi((ff + k))(ptq t ) + (Piqi){u(pi)u(q) + k) +pi(q~f + k)) 

-Pt(u(qf + k)uj{ql) + (q~) + k)qf) (uj(qf)cj(pi) ~ q"fpt) 

5(cj(q}) + u(qf + k) + u)(k) - 2uj(pl)) 



(q f - Pl ) 2 2q f (h - k){q f - Pl ) 2 2q f {l 2 - k) 



(5.32) 

Now we are left with the integration w.r.t. qj. In order to be able to do this we 
have to rewrite the argument of the delta distribution showing the dependence on 
explicitely. We want to use the formula 



n 



(5.33) 



where the simple zeros of the function /. For that purpose we define 



/(qf) = yqf + ™ 2 + y(q} + k) 2 + m 2 + - 2cj($ 



Vqf 2 + m 2 + Vqf 2 + k 2 + 2q f k cos 9 + m 2 + k - 2uj(pl) 



(5.34) 



where 9 = <(q),k) and we have used the notation k = qf = \q~f\. In order to 
simplify the notation in the following formulas we write u> = 2u>(pi). The zeros of 
the function / are given by 

q f U(k,9) = A(k,e)+B(k,e) 

q f (2) (M) = A(k,9)-B(k,9) 

where we have introduced 

.(2k-.)kcos0 



2(uj - k) 2 - 2k 2 cos 2 9 



(ui - k)^/uj A - 4w 3 k + 8wkra 2 - 2k 2 m 2 + 4cj 2 (k 2 - m 2 ) + 2k 2 m 2 cos 29 
S( ' ' = 2(w - k) 2 - 2k 2 cos 2 9 

(5.36) 

Clearly these solutions are itself functio ns of the momentum k and the angle 9. If 
we consider the limit of these functions ( 5.35| ) as k goes to zero we obtain 



limq f « =± Pi , j = 1,2 



(5.37) 
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where pi = \pl\. The derivative of / with respect to qf is given by 
± = qf + qi + kcos^ 

aqf ^/m 2 + q f 2 ^/k 2 + 77i 2 + q f 2 + 2kq f cos 

By inserting the two zeros of / into (5.3S) we observe that, by taking the limit 
k — > 0, we get a finite result, namely 

dfM 



lim 



qf=qf 



k^o dqf 

Now we can rewrite the cross section as 

1 (2tt)- 5 k 6 
a — — - 



= ±^,i = l,2 



(5.39) 



4 /((k%) 2 -™ 4 ^ 



dH 1 d%g(h)g(l 2 y 



d 3 k qf 2 dqfdQ 



2oj(k) 2w(g/) 2u%/ + fc) 



(uj(qy)uj(q) + k) + q f (q) + fc))Oi?i) + (uj(q})uj(ql) - q}ql)(uj(pl)uj(q} + k) 

+ Pi((ff + k)) - (u}(qf)w(fi) - q)pl) (u)(q f + k)uj(ql) + (q) + k)qt) 

- (uj(pl)uj(q f + k) +Pi(q~} + fc))(Pift) _ (PiQi)(w(pi)w(qf + k) +fi{q~f + k)) 

+ pf (w(q~f + k)u}(qi) + (q~f + k)qt) (u(q f ) 2 - q f 2 ) - (u(q))uj(q f + k) 

+ Qf(0 + k )) (PiQi) - - QjQi) (u{Pi)u{ff + k)+ pi(0 + k)) 

+ M<7/Mk) - <ffP~i) (M<7/ + k)u(qi) + {qf + k)qt) + (uj(pl)uj(q f + k) 
+ Pi(q~f + k))(p t q t ) + (Piqi)(u(pi)u(q) + k) +p l {q} + k)) 



-i 2 



°i (">(?/ + k)oj{ql) + (q) + k)qt) (o>(?/) w (Pi) - qjfi) 



(q f -pi) 2 2q f (h ~ k){q f - Pl f2q f {l 2 - k) ^ |/'(q f )| qf=qf «, | ^ ^ 

(5.40) 

Before we do the integration w.r.t. qf let us denote by i*i(qf , k, kcos0) the expres- 
sion which comes from the tensor part in the numerator, i.e. 



Fi(qf, k, kcos( 



M<?/M<7/ + k)+ q~}(q~f + fc))(pi&) + (uj(q})u;(cil) - q)ql) 

(uj(fi)cj(q~) + k)+ pi(q" f + k)) - (u}(q"f)ui(pi) - q}pl) (aj(q f + k) 
u{ql) + {q~f + k)ql) - (uj(pi)u){qf + k) + Pi(q~f + k))(p l q l ) - (p l q l ) 
(v(pi)u}(q~f + k) +Pi(q~f + k)) + p 1 2 (uj(q'f + k)u(ql) + [q) + k)ql) 
(uj(q}f - q} 2 ) - (uj(q})uj(q}- + k) + qjiqj + k)) (p^) - (uj(q}) 

~ q^cii)(v(Pi)v(q~f + k) +p~i(q~} + k)) + (uj(q~f)u(pi) - q}pl) 
(u>(q~f + k)Lu(ql) + (q} + k)ql) + (u;(jh)uj(0 + k) + Pi{q~f + k)) 
(PiQi) + (Pift)(w(pl)w((7} + k) +Pi(q~f + k)) - pj(uj(q~f + k) 



-i 2 



+ i<ff + k )<ji) (^(q})^{pi) - <TfPi) 



(5.41) 
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Then the expression for the cross section is given by 



a= 4 



1 f d % d %g( h m,r f d3k qf2dqfdn 1 



^((K%) 2 - ™ 4 ) J ' - - J 2u,(k) MQf) MQf + k) 

Fi(q f ,k,kcos0) V 1 Af _ {j) 

(5.42) 

From this expression we can identify the differential cross section 4°.. We do the 
qf-integration and obtain 

da 1 (2^)- 5 k 6 f . . w , f d 3 k 



(q f (i) ) 2 Fi(q f W),k,kcosfl) 



2y / ((qf°' ) ) 2 + ™ 2 ) 2y / ((q f 0)) 2 + k 2 + kq f W) cos<9 + m 2 ) 

/ 1 V (5-43) 

\m 2 - 2 v /(q f W) 2 + m 2 uj(pl) + 2q f W) p; cos a + p 2 / 

1 

2[\/(qf^) 2 + m 2 Z? - q f O')li cos/3 - A/(q f «) 2 + m 2 k + q f «kcosfl] 

1 

2[V(qf (3) ) 2 + m 2 ?5-q f wi2 cos 7 — 

V(qf W) ) 2 +m 2 k + q f 0')k cosy 

where we have introduced 

a = <(<?/, k) 
/?= <(#,£) 

(5.44) 

7 = <(<ff,h) 
li = i = l,2 

In the next step we omit all small quantities in the numerator, and we set k = 
in the denominators if there occur no singularities. Then the expression can be 
simplified to 

dn 4 jy—p m4) 2pi y y 2c( fc ) 

S 4^kF Fl(qf ° )|k = '°' 0) (m 2 - 2c(pl) 2 + 2q f O-)| k=oPiC o Sa + p?) 

1 



2[w(p5)ii - q f W|k=olicos/3 - w(pi)k + q f O)| k=o kcos0] 

1 

2[w(p 4 )/§ - qf (i) |k=ol2 cos 7 - ui(pl)k + q f 0')|k=okcos6>] 

(5.45) 
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Now one observes that the function F\ is independent of cos 8 and depends only on 
the initial momenta pi and qi so we can take it outside the integral. Then we get 

^> 2 (q f G) | k =o) / d%d%g(h)g(l 2 y 



da _ 1 (2it)- 5 k 6 
dtt ~ 4 



{{ViQiY ~ m 4 ) j=i 
d 3 k 



1 



(5.46) 



2oj{k) 2[w(pi)(/? - k) - q f O')| k =o(licos/3 - kcostf)] 
1 

2[co(pl){q - k) - qf0')| k= o(l 2 cos7- kcos^)] 
where the function F% is given by 



F2(q f (j) |k=o) 



^i(qf O) |k=o,0,0)| 



1 



m 2 - 2w(pj) 2 + 2q f O')| k= oPi cos a + p 2 y 

(5.47) 



If we use (5.37) we can write the differential cross section in the form 



5,^6 



da _ 1 (2tt)- 5 k 



dQ 16 '((Pift) 2 -"^ 



d 4 hd 4 l 2 g(h)g(l 2 y 



F 2 (+ Pi ) 



d 3 fc 
2w(fc) 



[w(pi)(Z° — k) — Pi Ci cos /3 — kcos#)] [u>(^)(Z§ — k) — p;(l 2 cos 7 — kcos6*)] 
d 3 k 



+ -FM-P0 / 



2w(fc) 



[w(pi)(/? - k) + pi(lx cos/3- kcos6<)] [w(pi)(Z§ - k) + Pi(l 2 cos7- kcosfl)] 

(5.48) 

In order to do the integration w.r.t. k we use an identity due to Feynman 



1 



dx 



a, 6 G 



a6 7 [ aa: + ^(1 — a: )] 2 ' 
We consider the first term in ( 5.48| ). Then we define the functions a and 6 by 

a= [uj(pi)(li — k) — pi ( 1 1 cos/3 — kcos#)] 

6= [w(pi)(ia — k) — pi (I2 cos 7 — kcos6>)] 

The integral becomes 

dx 



(5.49) 



(5.50) 



2w(fc) ah 2 • 

_ 2tt 
~ ~2 



[ojc + 6(1 - a;)] 

kc/k I dx f 

J-i [oas + 6(1 — a:) 



(5.51) 



We insert a and 6 and obtain 
7T I kdk J dx J 



dcos 6 



[— ui(pi)k + p;kcos# + (uj(pi)l1 — pJi cos (3)x + {w{pi)l2 — Pi I2 cos7)(l — a;)] " 



(5.52) 
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We set cos 9 = z and we observe that the integral w.r.t. z is of the form 

' _1 dz 2 



1 (cz + d) 2 d 2 - c 2 
sec where we've set 
c = Pik 

d= - u(pi)k + (u(pi)l° - Pili cosP)x + (u(pi)l° - pil 2 cos7)(l - x) 



Then we obtain for (5.52) 



= 2tt / kdk / dx 



P?k 2 + A(h, h, x)k + B(h, l 2 , x) 



Then this integral w.r.t. k is of the form 
ydy 



1 / 2 \ c 2 

ln(ciy + c 2 y + c 3 j 



dy 



ciy 2 + c 2 y + c 3 2ci 
where the constants Cj, £ = 1, 2, 3 are given by 

2 2 

ci = p 4 = m 
c 2 = A(lx,l 2 ,x) 
c 3 = B(h,l 2 ,x) 
In our case the discriminant A := 4ciC3 — c 2 is given by 



2ci J ciy 2 + c 2 y + c 3 



4tf 



- pik cos 0)x + (uj(pi)l 2 - Pil 2 cos7)(l - a;) 



= -4pf 



(w(pi)Z? - Pili cos/3)a; + (cj^)^ - pil 2 cos7)(l - x) 



so we see that A < 0. Then the integral (5.57) is given by 



ydy 



1 



ciy 2 + c 2 y + c 3 2ci 



ln(ciy 2 + c 2 y + c 3 ) 



f'2 



2ciV Z A 



In 



2ciy + c 2 



(5.53) 



(5.54) 



(5.55) 



where A(l\ 1 l 2l x) and B(l\,l 2 ,x) are given by 

A(/i,Z 2 ,a:) = - 2u(pi) (u)(pl)l° - Pi li cos/3)x + (w^JZj - Pil 2 cos7)(l - x) 

r n 1 2 

B(h,l 2 ,x) = ((J(pi)l° - Pilieos/3)x+ (uj(pl)l 2 - Pil 2 cos7)(l - x) 



(5.56) 
(5.57) 

(5.58) 



4uj(pl) 2 (oi(pi)J?- pJi cos (3)x+ (uj(pl)l°- Pi \ 2 cos (I ~ x) (5.59) 



2ciy + c 2 + V-A 

(5.60) 

Here we consider gravitons up to an energy ujq which is much smaller than the energy 
of the incident particles. Since every real detector has a finite energy resolution 
below which he cannot detect any particles we have to integrate over all these 
contributions up to the value u>q. Then the final result of the k- integration is given 
by the following expression 

kdk i r , 2 . 

ln(cicj + c 2 ujq + c 3 J 



1 



p 2 k 2 + Ak + B 2 Cl 



C2 



:hl 



2ciWo + c 2 



-A \2ciujq + c 2 



(5.61) 



ln(c 3 ) 



c 2 



In 



c 2 - x/^A 
^A 



C2 
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The second term in ( 5.4§| ) can be treated in the same way. T he o nly difference is 
the sign of pi in the definition of the functions a and b, see (5.5C), which has no 
consequences for the calculation of the integrals ( |5.53 ) and (5.57). 



5.3. Adiabatic limit. In the preceeding subsection we've found the differential 
cross section for bremsstrahlung in a scattering process of two massive scalar parti- 
cles. Now we want to discuss the adiabatic limit. For the differential cross section 
we've found 



da _ 1 (2tt)- 5 k 6 
dfi ~ 16 




d i l 1 d%g(l 1 )g(l2)*F2(+Pi) 



(5.62) 





C2 - ' 






( 


C 2 + - 







where c\ — m 2 and Cj = Cj(Zi, Za, a?), i — 2,3. This expression, as it stands, is 
well defined and finite due to the presence of the testfunctions. To obtain a phys- 
ically relevant result we have to remove this cutoff and therefore test the infrared 



behaviour of (5.62). We will show below that part of this expression becomes sin- 
gular in the adiabatic limit. In the adiabatic limit in momentum space we let the 
testfunctions tend to a delta distribution in it's argument. This will be done in the 



same way as in subsection 5.1, see (5.14). We scale the argument of the testfunc- 
tions with a parameter e in the integrals of ( 5.62| ) and consider the limit e — ► 0. 
The functions C3 and A are homogeneous of degree 2 in the variables l\ and l 2l i.e. 



c 3 (eh,sl 2 ,x) = e 2 c 3 (h,l 2 ,x) 
A(eh,el 2 ,x) = e 2 A(li,l 2 ,x) 



(5.63) 



while the function c 2 is homogeneous of degree 1, i.e. c 2 (el\, el 2l x ) — ec 2 (li,l 2 ,x). 
With this in mind let us now look at the various logarithms in J5.62| ) . The argument 
of the first logarithm contains the constant m 2 ojQ so it stays finite if e goes to zero. 
The constant in front of the second and the fourth logarithm has equal powers of 
e in the numerator as well as in the denominator so it tends to a constant. The 
argument of the second logarithm goes to one, so the logarithm itself goes to zero 
and the fourth goes to another constant. It remains to consider the third logarithm. 
This one becomes singular in the adiabatic limit. This can be seen as follows: 



da 7T (2tt)- 5 k 6 
dn ~ ~ 8 



{{pm 



d i l 1 d%g(l 1 )g(l 2 rF 2 (+p ] )^- I dx 



2ci J 



In 



(u>(p~i)l° - Pill COS 0)x + (uj(pi)l2 - Pi I2 COS7)(l - x) 



(5.64) 
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If we concentrate on the singular part only we obtain by inserting the scaled test- 
functions 

d%d%g{l 1 )g(l 2 )* \n\(u(pi)ll - pJi co&0)x + - Pil 2 cos7)(l - x)\ 

= J dHid A l 2 ge{h)ge(l2)* ln\{u(pl)l° - p;h cos/3)x + (u(pi)1% - pil 2 cos7)(l - x)\ 
^d 4 hd 3 l 2 go( — )9o( — y^\(^(Pi)li ~ pJi cos/3)a; + (w(p^2 ~ Pil2COS7) 



e° J s e 

X (l-3f)| 

d 4 kid 4 k2go(ki)go(k2)* bx| (uj(pi)ek® — pieki cos@)x+ (wijfijek® — p,ek2 C0S7) 

x (l-af)| 
= (27r) 4 ln|e|+o(l) 

(5.65) 

where we've introduced fcj = j. So we get a logarithmic divergence if e tends to 
zero. The origin of this singularity is the massless graviton propagator as can be 



clearly seen from the above calculation. Since all the other terms in (4.14) have the 
same structure of propagators we would get the same logarihmic divergence as in 
our example. Furthermore there can be no cancellation between these terms since 
they all have a positive sign. So we can omit the discussion of the other terms. 

This result is similar to the case of quantum electrodynamics || where a log- 
arithmic divergence in the differential cross section was obtained as well in the 
brcmsstrahlung contribution to the scattering of an electron due to an external 
classical source. As is well known from QED, the infrared divergence in the 
brcmsstrahlung must cancel against contributions from radiative corrections, oth- 
erwise scattering theory would break down. 

Appendix 

In this appendix we give the definition of the singular order to of a numerical 
distribution t S <S(R 4 ). The singular order measures the singularity of the distri- 
bution t at the origin ||. Let t G 5(R 4 ) be given. Then the distribution t has 
scaling degree s at x = 0, if 

s = inf{s' € R\X S t(Xx) —5 in the sense of distributions} 

The singular order is then defined by u> := [s] — 4, where [s] is the greatest integer 
less or equal to s. It's an easy excercise to see that differentiation increases the 
singular order of a distribution, i.e. 

co(D a t) =io{t) + \a\ 

where D a = d^l ■ ■ ■ d^i for every multi-index a = (ai, . . . 0,4) with \a\ = a% + . . . + 04. 
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